LONG TIME EXISTENCE OF SMOOTH SOLUTION FOR THE POROUS MEDIUM EQUATION IN 
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i Abstract. In this paper, we are going to show the long time existence of the smooth solution for the porous 

medium equations in a smooth bounded domain: 

■ (ttf = An'" in Q x [0, oo) 

(0.1) I u(x, 0) = u Q > in Q 

\u(x,t) = forxedQ 

<q 

where /;; > 1 is the permeability. The proof is based on the short time existence of C?" -smooth solution, the 
global Q -estimate, the Holder estimate of divergence type degenerate equation with measurable coefficients and 
C]' 1 ' -estimate of mixed type equation with Lipschitz coefficients. 
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1. INTRODUCTION 



' We consider in this paper the initial value problem for the Porous Medium Equation(PME) 

> ■ 

\Q ■ f u t = Au m in Q 

(1.1) L(x,0) = u >0 inQ 
\^ \ \u{x, t) = for x 6 dQ. 

posed in a bounded domain Q with the range of exponents m > 1, with initial data Uq nonnegative, integrable 
and compactly supported. 

In 1 6], Kim and Lee dealt with the short time existence of solution to PME in a bounded domain. More 
precisely, the main result in their paper is that, if u is a solution to ( jl.H and / = u'", then, in some regularity 
conditions on the initial data f° = u™ and its first and second derivatives, the solution of the degenerate 
equation 

(1.2) fi = mf^r A / = mf a/, if = u m , a = 1 - 

exists on a short time interval [0, T) and / e cf } (O) on [0, T), i.e., the solution and its first and second 
derivative are Holder continuous with respect to a suitable Riemannian metric s under which the distance 
between two points X\ and X2 in Q is equivalent to the function 

|xi - x 2 \ 



d( Xl ) a + d(x 2 ) a + Lti l*i - *2\ a 

with d = d(x) denoting the distance to the boundary of Q. 

In this work, we will show that, under the same assumptions, the solution f of (11.21 exists on the long 

time [0, oo) and it is also of class (ff , i.e., assuming that the initial value f° is strictly positive in the interior 
of a domain Q C R", with f° - on the boundary, and denoting by d the distance to the boundary of Q, we 
will obtain the following result. 

Theorem 1.1. If the functions f°, Df° and d"D 2 f°, restricted to the compact domain Q, extended continuously up 
to the boundary ofQ, with extensions which are Holder continuous on Q of class O'(0,),for some y > and Dfo + 

l 
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along dd, then the initial value problem 

'f t = ±f«Af, (x,t)eOx(0,oo), 
(1.3) \f(x,0)=f(x), xeQ 

/(x, f) = 0, xedClx [0, oo) 

admits a solution f which is Cf y -smooth up to the boundary, when < t < oo. 

As in the Section 4 in |6|, the coordinate change, (z = f(x', x„, t) — > x„ = /z(x', z, f)), converts domain 

DeR" => D 6 R"" 1 x M + 

and the equation <1.2t to 

1 + IV^P 



(1.4) h 



A x ,h + 



h z 
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In addition, Ti^., (z = 1, • • ■ , n - 1) satisfies 

(1.5) w t = z a V k (a kl Vnv) , (k, I = 1, ■ ■ ■ ,n) 
where 

a nn = 1 + lV 7 x ' hl \ a k '" = -^, a" k '=0 and a k ' v = 5 kl , (k'J = 1, 1). 
hj K 

Since the solution f of (|1.31 is strongly related to the solution h of {1A) , by the Theorem 1.1 in (6), the 
solution exists on a short time interval. Let (0, T) be the maximal interval of existence for C^-solution h 
and suppose that T < oo. Then, solution does not belong to the space at time t = T anymore. However, 
the existence theory gave h e W 1 - 2 . This implies that the coefficients of the equation dl.4t are only measurable 
and bounded. Thus, if it could be shown that h 6 c]' r , then the coefficients in lll.4|l would belong to . 
Then, Schauder estimate in |6| provides Cf y regularity on h at t = T and we can get a extended interval 

[0, T'), (T < T') in which / e C^' 7 . This contradicts the maximality of T. Therefore, T must be infinity and 
the Theorem ll.ll follows. Hence, the missing step for the regularity problem to be solved is 

(1.6) h e W 1 ' 2 ^/ie C 1 /, (0 < y < 1). 

In this paper, we show the long time existence of solution f - u m by solving the missing step l ll.6|l 
The paper is divided into five parts: In Part 1 (Section 2) we review the metric ds which controls the 
diffusion and state theorem for the short time existence of / = u m in |6). In Part 2 (Section 3) we deal with 
global estimates, including Gradient estimate and Non-degeneracy. In part 3 and 4 (Section 4 and 5) we 
establish the Holder estimate for the solution of degenerated parabolic equation with divergence type and 

Cg' estimate for the solution of mixed equation with Lipschitz coefficients. Finally, in Section 6, we prove 

the long time existence of the solution / of dl.21 which is in Cg '' (□). 

2. Preliminaries 

The diffusion in \1A\ is governed by the Riemannian metric ds where 

dx? H + dxl 

ds 2 = — i -. 

The distance between two points x l = (xj, • • • ,x\) and x 2 = (x 2 , • • • ,x 2 ) in this metric is a function sLx^x 2 ], 
which is equivalent to the function 



■ r 1 2l ^ 1=1 

\x\f-\ - 



X 2 \ 



\x\\^+K\i+n:l \x)-x]\i 
in the sense that 

s <Cs and s < Cs 
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for some constant C > 0. For the parabolic problem we use the parabolic distance 

s [(x\ h), (x 2 , t 2 )] = s [x\ x 2 ] + Vlfi - hi 

In terms of this distance, we can define Holder semi-norm and norm of continuous function g on a compact 
subset A of the half-space {(x\, ■■■ , x n , t) : x n > 0}: 

„ „ g(Pi)-g(Pi) 
Wtf W = su P ^vT' 

P l5 tP 2 €A S{t / \-t J 2) r 
I^CfCA) = l^llc°(A)+ll^ll H r (A) , llgllc°(A) = Slip P€A \g(P)\ . 

With these norms, the space cf y (A) is the Banach space of all such functions with norm: 

n n 

IIS"cf (A) = ^Cl(A) + E HfeUcftA) + H^llcffA) + E IK^/HcTCA)- 
!=1 l<i<j<n 

Imitating the case where the operators are defined on the half-space {(xi, • • • ,x n , t) : x„ > 0} we can define 
the distance function s in O. In the interior of Q the distance will be equivalent to the standard Euclidean 
distance, while around any point Xq e dCl, s is defined as the pull back of the distance on the half space 
{{x\, ■ ■ ■ ,x n ,t) : x„ > 0) through a map (p : {(x\, ■ ■ ■ , x n , t) : x„ > 0} — > O that straightens the boundary of Q 
near x . 

It can be easily shown that the distance between two points Pj and P 2 in Q is equivanlent to the function 

- m p \ IP1-P2I 

s(Pi,P2) = ; 

d(p 1 r + d(p 2 r + Li:i\Pi-P2\ a 

with d = d(P) denoting the distance to the boundary of Q. The parabolic distance in the metric is equivalent 
to the function 

s [(Pi, h), {Pi, fa)] = s [Pi, P 2 ] + Vlfa - fal- 

Suppose that A is a subset of Ox [0, 00). As above, we denote by C7 S (A) the space of Holder continuous 

functions on A with respect to the metric s and by cf~(A) the space of all functions w on A such that u>t, 
Wj and d a w; j , with i, j e {1, • • • , n} and with d denoting the distance function to the boundary of Q, extend 

continuously up to the boundary of A and the extensions are Holder continuous on A of class Cl (A). Then, 
they are both Banach spaces under the norm IMI r ~ (A1 an d 



\\W\Uy = \m, c 



f(A) + E IKIIcfW + l|lytll Cr(A) + E H^^/UcffA)- 
(=1 l<i<j<n 



Under this metric, it is known that the problem l|1.3t has the cf } solution for a short time interval (0, T). we 
now state the short time existence for ([1.31 1. 

Theorem 2.1 (Theorem 1.1 in J6]). Let d be the distance to the boundary Q. If the functions f°, Df° and d a D 2 f°, 
restricted to the compact domain Q, extended continuously up to the boundary o/Q, with extensions which are Holder 
continuous on Q of class O' (CI), for some y > and Dfo + along dCl, then there exists a number T > Ofor which 
the initial value problem 

C f = mf a Af, (x, t) e Q X (0, 00), 
\f(x,0) = f°(x), xeCl 
[f(x, f) = 0, xedClx [0, 00) 

admits a solution f which is C^' 7 -smooth up to the boundary, when < t < T. 
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3. Global Cj| -estimate 

This section is devoted to prove some properties of the solution / to jl.3i , including global estimates, 
gradient estimate on the boundary, etc... . To get them we construct sub and super solutions to the problem 
111. 31 . Now, we first deal with the L°° estimate of /. 

Lemma 3.1 (L°° -estimate). Let d be the distance to the boundary ofdD, and let f be the solution of il.3i with initial 
data f° such that 

f°, Df°, d"D 2 f° : continuous up to the boundary dQ. 
There are constants c > and C < oo such that 
(V 



ll/llL°°(n) 

(2) There is a ball e Q such that 



- min {(n^' ll/ollL " (Q) }' ( a = 1 ~^) 



mif(x,t) > 



A > 



(1 + f) 1/£l ' 

Proof, i) We take a ball Br = Br(0) of radius R strictly containing Q, especially we assume that QcBs, and 
consider the function z(x, t) defined in Br X (0, oo) by 

(l + t)= 

for suitable constant A to be chosen presently. Since z is positive in Br x (0, oo), we have 

f(x,t) = 0< z(x, t) on <9Q x (0, oo). 

In addition, if we choose A larger than R2 °° <Q) , we get 

fo(x) < z(x,0) in O. 
Finally, we will obtain the inequality mz a hz - z t < in Q whenever 

R 2(i-a) - 1 
2mna 

With this choice, the comparison principle implies that 

f{x, t) < z{x, t) = A(R2 ~ W2) < = — ^ in Q x (0, oo). 

+ (i + f)« (i + 0« 

On the other hand, one can easily check f(x, t) < ||/oIIl~(0) by comparison principle. Hence the proof of (i) is 
finished. 

ii) We first suppose that the initial data Uq = fj satisfies the following condition, 

(m - l)Au™ < -uq Vx e O. 
Then, the initial data uq is controlled from below by the solution g(x) of 

' (m - l)Ag"'+g = in O 
g =0 on dQ, 

i.e., 

U Q (X) > g(x). 

In addition, the function 

(1 + 0^ 



LONG TIME EXISTENCE OF SMOOTH SOLUTION FOR THE POROUS MEDIUM EQUATION IN A BOUNDED DOMAIN 5 

is also a solution of porous medium equation with the initial data g. Hence, by the comparison 
principle, we have 

u{x,t)> * K \ . 

(l + t)& 

Therefore, for any ball B§ e Q, we have 

M xeB g(x) 
inf f(x,t) > —z—. 

.T€B S / V (1 + t)* 

Next, we denote by Q„ the set 

Q„ = {x e Q : (m - 1)am™ > -« | 
and we assume that O„ + 0. Let's define the function u Q such that 

«o = Mo on O„ 

and 

(m - 1)am " > -Mo on Q. 

Then, by comparison principle for elliptic equation, «o < uq in O. We also let Ti be the solution of ^1.1) with 
the initial data Uq being replaced by Uq. Then, by comparison principle for parabolic equation, we have 

u > u in Q. 

On the other hand, by the Problem 8.1(i) in (9|, u satisfies 

_ u 

Ut > . 

f - (m-l)(l + f) 

Thus, by the Gronwell's inequality, we have 

-I A 

M X,f > — • 

(1 + f)* 

Since Q„ is open set, there exists a ball B 6o e Q„ . Hence 

— , - Mn(x) Mo(x) 

m x, f > m x, t) > = on B 6o . 

Therefore 

inf YSB u"' (x) inf xeB f (x) 
inf f(x,t) = inf u m (x,t) > ^4 — = ^4— 

«EB 8 / V x<EB So (1 + 0* (1 + 0* 



, which implies the conclusion. 



Next property is the gradient estimate which will play an important role to show C s ' y continuity of 
solution. 

Lemma 3.2 (Gradient estimate). Let d be the distance to the boundary of 3D, and let f be the solution of jl.3i with 
initial data f° such that 

f°, Df°, d a D 2 f° : continuous up to the boundary dQ.. 
There are uniform constant < C < oo such that 

ll/ollcMQ) 



(3.1) IIV/IHq) < Co 



(i + ty a 
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Proof. By the Theorem 12. II there exist constants to > and Ci > such that 

(3.2) |V/(z, 01 < Q VieQ, < f < t - 

Thus, ip.ll holds for a short time. 
Let ip be the solution of 

1 i/'(x) = on <?Q. 

Since < |V/o| < °o on <9Q, we can select constants 0<Ci<l<C2<°o such that 

cnp(x) < fo(x) < c 2 ip{x), in O. 
Then, by the comparison principle, we get 

<f(x,t)< f, xeQ. 



(l + f)« (1 + 0" 

Since Hi/'Mllcvo) < 00 / we also have 

IV7 r/ , \ i ^ Pzll^Wlb(aq) „ c 2 ||/o(x)|| c i (Q) 

(3.3) Vf(x,0 < r — - ^ r^- < °°/ V0 < f < oo, x e ,9Q. 

(l + t) = d(l + f)« 

Denoting by Q CT , for a > 0, the set 

Q a = {x e Q : dist(x,<?Q) > a). 
Then, by | |3.3l l, there exist constants oq > and C3 > such that 

(3.4) |V/(x,0l < / < 00 VO<t<oo,xeO\Q (Jo . 

(1 + 0° 

To finish the proof, let us define the family of rescaled functions 

f k (x,t) = k^f(x,kt) 

in Q„ with parameter k > 0. Then, they are again solutions of 

(3.5) ft=mf a Af. 

Since / > in Q O0 , the coefficient mf" is bounded above and below. Hence the equation d3.5l l becomes 
uniformly parabolic and the solutions are smooth with derivatives locally bounded in terms of the bounds 
for /. We conclude that there exist a constant C\ > such that 

\Vf k (x, 1)| < c 4 xeQ ao . 

This means that 

(3.6) \Vf(x,k)\<^T xeCl O0 . 

k« 

For some constant C5 > 0, putting t = k and C4 = Csll/ollc^o) m inequality d3.6b . Then 

fr, rn rnr/ x\i c sll/olb(Q) „ 

(3.7) \Vf(x,t)\ < ; «Q„, 

in 

Hence, we get 

(3.8) |V/(x, f)l < C6 " / ° llcl( i n) , V* > t x e Q lJo 

(1 + £)» 

where c 6 = c 5 (l + . By (12), (ED and 0, (E) holds for all f > 0. □ 

Finally in this section, we will show the non degeneracy of solution f to 111. 31 near the boundary. This 
estimate quarantees the uniformly ellipticity of coefficients a kl in ( |1.5t . 
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Lemma 3.3 (Non-degeneracy estimate). Let d be the distance to the boundary of dQ. and let f be the solution of 
Jl .3I > with initial data f° such that 

f°, Df°, d a D 2 f a : continuous up to the boundary dQ. 

There are uniform constant < C\ < oo such that 

(i +° t )i/ a < H v /Hl~(5o) 

Proof. By the Lemma l3~Tl (ii), there exists a ball B,5 c O such that 

inf f(x,t)>— ^-p 

* eB * Q (l + f)« 

For fixed t > 0, let v be the solution of the problem 

hv(x, = in Q\B 6o 

v(x, t) = c on dB 5o 

v(x, t) = on dQ. 

Then, the function V(x, t) = D ^ x,t \ satisfies 

(l+t)* 



In addition, we have 



and 



V a AV -V t = - > 0. 

a(l + f) 1+ « 

V(x, t) = = < f(x, t) on dB 6 

(l + t)i (l + 1)« 

V(x, t) = V ^ X,t \ = = fix, t) on dQ. 
(1 + t)- 



By the comparison principle, we have 

f{x,t)>V{x,t) inO\B flo . 

By the Hopf 's inequality for the harmonic equation, there exists some constant C\ > such that 

dv „ ^ 

— < -Ci < on dQ. 
dv 

for the outer normal direction v to d£l. Therefore, 

\\Vfh<°(da) > \\W\\ L ~ m = M > ^~ 



(1 + f)« (1 + f)« 

and lemma follows. □ 

4. Holder estimate I 

In the previous section, we obtained the global C.J regularity of solution / to the problem ( 11. 31 . By the 
relation between / and h, we can say the same story on the solution h of (1.4) . Hence, we can have basic 
informations for coefficients a kl and solution w in \1.5\ , including boundedness of solution and coefficients. 
With this basic properties, we devote this section and next one for solving the missing step, (1.6) , for the 
regularity theory. 

Let H be the half space {x = {x\,- ■ ■ ,x n ) e R" : x„ > Oj. We are going to show Holder estimate on a 
solution zv of the equation 

(4.1) w t = xy t UNjw) + x a n g in H 

for a forcing term g. Assume that the coefficients a ll (x, t) are measurable functions and satisfy 

(4.2) A|£| 2 < a'i{x, t)^ < A|£ | 2 , (i, j = 1, ■ ■ ■ ,n). 
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In addition, we suppose that the forcing term g satisfies 
(4.3) \g(x)\ < C\w(x)\ 

for some constant C < oo. 

For the Holder estimates of the solution w to 1)4. 1) , we need the following two inequalities. The first one 
is a weighted version of Sobolev's inequality. Let Cj^(H) be the space of restriction of functions in C^°(R") to 
H. 

Lemma 4.1 (See Theorem 4.2.2 in (8)). Let 1 <p < q < oo, s > — jj and we assume that a < 1 satisfies 

1 1 -a _ 1 
q n p' 



Then 



J x%u p dxj <c|j" cbcj" . 



for the closure of C™(H). 

The next inequality is the parabolic version of Sobolev's inequality with weight: 

Lemma 4.2. Let Q = B x (a, b) is cylinder in H x (0, 1) and f e C™(Q). Then there are constant C > such that if 
I _ "+2-2« we j iave 

1 n—a 

r b r \( r \ r b r 

x~ n a \f\ 2h dxdt) <C sup x~ a \f\ 2 dx + \Vf\ 2 dxdt 

J a JB I \ (a,b) Jb ) J a JB 

Proof. By the Holder inequality 



x- a \frdx\ <Q \ x- n a \f\ z dx\ \\x- a \f\—.dx 



for some costant C\ > 0. By the Lemma [4. 11 it follows that 



f^ff* dxj <C 2 |£x-"|/| 2 dxJ k ^IV/pdxJ' 1 . 



for some costant C2 > 0. Now, taking the the l\ power and integrating in {a, b), we get 



f r x- a 

Ja JB 



\f\^dxdt\ <C 3 



fsup f x-"\f\ 2 dx + f f |V/| 2 rfxdf 



for some constant C3 > and the lemma follows. □ 

Define the balls B r and of radius r around x = Xq to be the sets 

B r (x ) = {x e R" : \x - x \ < r] and B^(x ) = B r (x ) n jx„ > 0|. 

We let B r be the ball around the point x = 0. We define the round cubes Q, of radius r around (x, f) = (0, 1) 
to be the sets 

Q r = B r x (1 - r 2 ~",l). 

We also define the general round cubes: 

Q r (x, t) = Q r + (x, t) and Q r + (x, t) = Q r (x, t) n {x„ > 0}. 
Let us give the Harnack inequality first. 
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Lemma 4.3 (Harnack's Inequality). Let w be a solution of equation (|4,l|l defined in B| X (0, 1) with conditions l|4.2Jl 
and g3J. Let \w\ < Mon Q\ andh = w + M + 1 > 1. Then, 

max h < Cmhxh 

Q+(o,-§) Q\ 

5 8 

Proof. Let <3 + = S + x (S2, sj e H x [0, oo) be an round cube in Q| and we take <p(x, t) = if(x, t)h 5 (x, t) as test 
function, where rj e C°°(<2 + ) with r\ = on d p Q + n jx„ > 0}. Since zt> is solution of l(4j) in Q+, /z satisfies 

(4.4) h. = y i {a i 'V j h)+g inH. 

Hence, multiplying l|4.4t by the test function <p(x, t), we arrive at 



^ 8A(1 

|l + or jjq 



+2 rr rj 2 ^ 6 igi ^+ rr x^p^^dxdt 

J J<0+ |1 + 0| JJ(0+ 



;q+ -1- jjq 

for S>0orS<— 1. Use of Young's inequality yields 



(4.5) 



sup f x-"i 1 2 h 1+6 (-,t)dx+ ff kUh^f dxdt 

i<f<s 2 J'B JJQ+ ' ' 

~ Ci IT ( |V7?|2 + x "" r? ' 7?f| ) (^) 2 dxdt + ff i]2 ^ 6 ^ dxdt 



for some constant C\ > 0. 

i 

We first consider the case 6 > 0. Let t\ and r 2 be such that \ < r 2 < t\ < (j) 2 ' , choosing rj in such a 
way that r]{x, f) = 1 in Q+ (0, -§), »j(x, = on d p Q+ (0, -|) n H, < tj < 1 in Q+ (o, -|), |Vf/| < ^- and 
|Tjf| < (ri _^ )2 _„ for some constant C 2 > 0. Then, for <3 + = Q + _j_ (o, -§) = B + x (|, |), 

(4.6) ff Q {m 2 + P) 2 < - _ c ; z|2 _ it Jfj^ *,r &o. 

for some constant C3 > 0. On the other hand, by ( 14.3b , the second term of the right hand side in ((43) is 
changed to 

(4 7) f (rr 2 /; 6 ) M dxdt < C 4 f (ifh 1+6 ) \w\ dxdt < C 4 f h M dxdt 

for some constants C4 > 0. Applying J4.6I I, J4.7II to 1 14.51 1, we can get 



sup f x-"ifh 1+6 (-,t)dx+ ff \y{qh^)f dxdt 

L<rf<5 J'B+ JJQ+ ' ' 

- ff x- a (h^f dxdt 



i<f<5 Ja 

(4.8) 

L 5 



In - r 2 | 2_ 
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for some constant C5 > 0. By the Lemma [4.21 there exists h > 1 such that 

1 

ff x-"(h^f h dxdt) 



(4.9) 



sup r x - a h i+6 (;t)dx+ rr ivrwi' 

i< f <5 Js+ JJq+ 1 1 



for some constant > 0. Then, by J4.8II and 114.91 , 



(4.10) 



rr & mt d X dtf < , C7 l2 r (m a 



dxdt 



for some constant C7 > 0. Taking the irjth root on each side of J4.10L we have 

C 8 

PII L a + ,, (l(Q . ( ,_ i)A7 .) < - _ IWIl-(q+ (o,-§K») 

for some constant Cs > 0. Define yi = po > 0, y, = Ziy,-i and r, ■ = g + ^(|) 2 " - |j • Then, the previous 
inequality becomes 



2-a 

."1 



"^(^.("-IM - C8 (q) ' t 8 ')''' 



Therefore, iteration yields 
(4.11) 



max h < G)\\h\\ 



, (C 9 > 0). 



For the other case 5 < — 1, we also take (p(x, t) = rj (x, t)h (x, f) as test function, choosing 77 in such a way that 
rj(x, t) = 1 in Q+, jj(x, = on <9 p Q+ n H, < r? < 1 in Q+, |Vrj| < ^ and |r] f | < for some constant 

1 

C10 > and | < r2 < ri < (|) 2 • Then, by the similar computation above, we have 

Cu 



PIIl< 1+5 >'o(Q+ ,V) - 



In - r 2 |i+« 



for some constant Cn > 0. Define yo = ~Po < 0, y* = hyt-i and r, ■ = | + 
inequality becomes 



1 \! 
1 j_ / / 1 ^ 2 -° 



I j . Then, the previous 



n 



Lw+i(cy. +l ,^«) > Cu (-] (8')''' ||/z||Ln ( Q+,. v >). 

, (Cu > 0). 



min/i > C12PH 

^1 IrPo 



Therefore, iteration yields 
(4.12) 

I (1)*= 

1 

To get the harnack's inequality for h, we finally examine the case 6 = -1. Let < p < (ly ~" . We take 
cp(x, f) = rf(x, t)h~\x, t) as test function, where 77 e q°(B+ X [0, 1]), i](x, f) = 1 in Q+ supp(r?) c Q + s c Q+ 

and Vj] < |)]t| < ^7 for some constant C13 > 0. Then, since if satisfies l l4.lt in , we can get 

( 



— r ff |Vlog/;| 2 dxdf 

{\x-"(B + p )\ JJq; 
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_ 2(n-a) 



for some constant Cu > 0. Let U = log/z and p = 2 . Then, by the Lemma [4.11 there exists constant C15, 
Ci6 > such that 



— i— , rr x-"\u-u Q; \d X dt<\ x -"(B;p- f (r x - a \u-u Q fd X 

X n a {B + p )\JjQ+ Jl-p2-«[JBi 



dt 



(4.13) 



< C 15 p 



¥ I if 



\VU\ 2 dx 



dt 



< Ci6p 2 2 - Ci6p 2 < Ci6- 

The parabolic version of John and Nirenberg Lemma for BMO (see [1 . 4]) yields that there exist two positive 
constants Pq and Cn such that 



(4.14) 



^ J x- a eP° u dxdt 



Then combining this with ||4.11| | and d4.12b , we get 



1 1 . x "' e 



-P° u dxdt 



< c 



17- 



max ft < Cis min h 

q\(o,-i) et 

H 8 

for some constant Cis > and lemma follows. □ 
We can now state the first main result of our paper. 

Theorem 4.4. Let w be a solution of equation defined in X (0, 1) with conditions d4.2b and <|4.3t . Suppose 
that 

max \iv\ < C\ w \ < oo. 

Qt 

Then, w is locally Holder continuous in X (0, 1) and 

IMb( Q +) < C (lMI L ~(Qn + l ) ■ ( s < !)• 

Proof. For h = zv + C\ w \, let m r = infQ+ h, M r = supQ+ h. Then h — m r + l,M r — h+l > 1 and satisfy the equation 
l|4.4)l in Q r . Applying the Harnack inequality (Lemma l4.3t to those equations, we get 

M r -Ml +1 = inf(M r - ft + 1) > - sup (M r -ft + l) = i(M r -mx +l) 

C Qt(o,-f) 



and 



Hence, 



/// < 



i 1 

- m r + 1 = inf(ft -m r + l)> — sup (ft - m r + 1) = — (Mi - m r + i) . 

Q1 C Q+ £ (o,-f) 



Mr-mr<( §-4) (M r - m r ) + X) 
8 \C + 1/ V C + l 
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Let osc(r) = M r - m r . Then 

(r\ (C-l\ /N 2(C-1) 

By an elementary iteration, we get 

osc(r) < Cr y [osc(l) + 1] 

for some constant C > depending on C, i.e., h is Holder at (x, t) = (0, 1). Since h = w + C\ w \, we also have 
Holder continuity of w at (x, t) = (0, 1). Finally, by translating and dilating, we have 

IMbton ^ c" (iNli-(Qf) + i) , (s < i) 

and lemma follows. □ 

5. Holder estimates II 

Through the previous section, we have established that the gradient of h in | |1.4t with respect to Xi, 
(i = 1, • • • , n — 1), is Holder continuous near the boundary (9Q. However, one can easily check that h z is not 
a solution of (1.5) . Hence, solving the missing step (1.61 is not completed at this moment. Thus, we will 
devote this section for the Holder regularity for h z . 

Let iv be a solution of the problem 

( Wi = x" {rfiwij} in H 

(5.1) I w(x,t) = on x„ = 

[ zv(x, 0) = Wo(x) in H 

with initial value wq{x) e C 0,1 (H). Assume that the coefficients a l, (x, t) are measurable functions and satisfy 

(5.2) A|5| 2 < a'i{x, t)^ < A\E\ 2 , (i, j = l,---,n) 

for some constants < A < A < oo. 

To get the Holder regularity for h z/ we will first construct an important, for our purpose, barrier function. 

Lemma 5.1. LetQ = B r + x(0, \)andP r {x') = (x',r) e HforO <r< landletK = {x e B+ : \x, - (Pl{0)).\ < f , t = 1,- 
Tften, f/tere exzsfs a solution f > and h'me T(r) > such that 

'x*(i0fij)-ft = O (x,t)eQ 
f = <9,Q 
|V/| > c(T) > 5Bl n {x n = 0], t> T. 



(5.3) 

Moreover, at t = 0, / satisfies 



supp{/(-, 0)| c K and < /(x, 0) < m l K < oo 



(5.4) 

for a constant mo > 0. 

Proof. In this proof, we will use a modification of the technique of (TJ to prove the lemma. For each x' e ]R" _1 
such that |x'| < ^, let Pi = Pi(x') and consider 



(x,t) = 



and g(x, t) = e g(x, t + tq) for To > we can choose later. By a direct computation, we can get 



g f (x,f) = e~ M, £(x,f + T ) 



-M- 



(}\X-Pr\ 



and 



K = e 



-Mt, 



(X, f + T ) 



2(f + T ) (t + To) 2 

4^ 2 (x I -P i(i )(x / -P i ,J 



(f + T ) 2 



(i + j and 1 < i, j < n) 
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and 

Hence 
if 

(5.5) 



gh =e- Mt g(x,t + T ) 



2 p 4f(xi-P iA f 



£ + T 



<(a ii gi,)-gt>e- Mt g(x,t + T ) 



(f + T ) 2 

2jSAr" 4jSr 2 



(1 < i < n). 



M- 



T 



> 



2BAr" ABr 2 
M> — + Ar- 



Now, we define the function g e by 

gc = max{0,g-e } 

for small eo > 0. Since g e takes maximum value of two subsolutions to j5.1ll at each point in H, it is also a 
subsolution of i5.1\ . We are going to choose proper constants |8, To and eo so that 

(5.6) supp{& o (-,0)}cK and < g Cq (P r , 0) < m 
and there is a constant T > satisfying 

(5.7) supp{g eo (-,r)} = Bj(P i ). 

For d5.6b , we take sufficiently small constant To and large one jS >> 1 so that 

1 



1 ^ 

(5.8) < e and eo < t <m + e . 

(47TT ) 5 (4TIT0) 1 

Then, 

g(; 0) < e on B r + \X and e < g(P r , 0) < m + e . 
Thus, the condition | |5.6l l holds for sufficiently large constant To and jS. 

For | |5.7| |, we will focus on supply (■, t)}. Since g eo (x, t) is radially symmetric about the point P>., for each 
t > 0, eo-level set of g e (-, t) is a ball centered at Pi . Let's denote by R(t ) the radius of the e -level set of g cg (-, 0- 
Then, one can easily check that there exists time t = f such that R(t) is increasing on (0, t') and decreasing 
on (f, oo). Hence, at time t' > 0, R(t) satisfies 

d{R(t)} 



dt 



= 



at t = t'. 



To show that there exists a time T > such that 

supply (vT)} = 6^), 
we will show that R(t') > i with the following equation 



(5.9) 



-Alt- 



'«o =e [4TT(t + T )]2. 



Differentiating on both sides of J5.9II with respect to t, we can get 

e n 



(5.10) 



M + 



2(f + T )) 



T 2 M 



att = t'. 



By l|5.5b , the last term of l|5.10t is bounded below by 4r 2 . Hence, we can get 

R(f') > 2r. 

This immediately implies that there exists a time < T < t' such that 



R(T) = 
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and 05.71 holds. Since Il5.7t is true, one can easily check that there exists a constant c = c(T) such that 

Wg eo (P ,T)\>c(T)>0. 



Next, we define function G by 



( 



G(x, t) - e\ 



4 \x-Pi 



|2\ 



1 - 



1 - 



4 X"P: 



|2\ 



-Mf- 



inBj(P f )x[T,l). 



(47T(f + T ))i 

for some constant e! > 0. Note that if we choose £\ sufficiently small, then 

G(x,T) < g e (x,T), inBr(Pj) 

since, for some constants c\ and C2, 

\Vge(; T)\ > ci > and |VG| < c 2 < °o on <?B|(Pi) 
On the other hand, by a direct computation, we can get 



G t (x, t) = G 



-M- 



2(f + T ) (f + To) 2 



and 



Gh = G 



i*- p v)( x i- p u) 



(t + to) 2 

r 2 (t + to) 



(i + j and 1 < i, j < n) 



and 



Then 



(5.11) 



Gh = G 



2j3 4£ 2 (x,-P y ) 2 



t + T (f + T ) 2 

M 2j3r" 



Wxi-PLj 2 8 
r 2 (t + T ) r 2 



(1 < z < n). 



x«(a'iG,j)-G t >G 



2 T + T 16(T + t ) 2 



M 
2 



1 



4 x - P. 



|2\ 



r 2 (l + T ) 



inB,(P f )x[T,l). Let 
(5.12) 



„ 64A(1 + t ) „ 256A 

|S > — — - and M > 



Ar 2 ~ 15r 2 ' 

Then, the second term on the right hand side of l|5.11l l is positive. Hence, 

xtfi'Gijj-Gt >0 mB^PrJxlT,!) 

if 

(ABAr 1 * 8Br 2 256A) 

5.13) M>max ^ + _^_,±^£H 

1 T t 2 15r 2 / 
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Now we let / be a solution of 

(*%(j'fi j )-ft = (x,t)eQ 
(5.14) f(x,t) = (x,t)ed,Q 

[f(x,0)=m l K x 6 

Then, for sufficiently small £\ > and M in J5.13b and jS satisfying l !5.8ll and 115.121 , we have 

f>g e inB r + x(0,T] 

and 

f>G inBiiPrfxlTA) 

by the comparison principle. Since x' e 1R" -1 is an arbitrary point with \x'\ < fg,by the definition of functions 
g eg and G, the solution f has nontrivial gradients on dB\_ D [x n = 0), i.e., the third inequality of l|5.3l l holds. 

To complete the proof of lemma, we are now going to show the existence of solution / to l|5.3b with initial 
condition l|5.4t . We begin by constructing a sequence of approximate domian so as to avoid the degeneracy 
of the equation. We may simply put 

t„ > - \ . 
n) 



B+ n = B + r r\{x n 



We now solve the problem 

' < {^(fnh) ~ (fn)t = {x, t) e B+ X (0, 1) 

(5.15) • /»(*,*) = (x,t)e^„x(0,l) 

f n (x,0)=m Q l K xeB+ n . 

The maximum principle, which holds for classical solutions, implies that 

< f n < m in B+„ x (0, 1). 
Moreover, again by the maximum principle 

/„</„+! inB+„x(0,l), Vn>4. 

Hence we can define the function 

f(x, t) = lim f n (x, t), in Q 
n — »oo 

as a monotone limit of bounded non-negative functions. On the other hand, we also have 

fn(x,0) < —^ x n inB ^H- 

Hence, < / < Mx„ for some constant M > and f is continuous up to the boundary. Since /,, is a classical 
solution of 1 15.151 1, / clearly satisfies d5 -3b and ( 15 -4b . Hence, we complete the proof. □ 

With this barrier function, we are going to show the second main result for our paper: C 1 '' regularity of 
solution w to H5.lt . 

Lemma 5.2. Let wbe a solution of 05. 11 1. Ifw is in C 0,1 (H), with respect to space variable x, then, it is also in C 1,y , 
(0 < y < 1) near the plane {x„ = 0], i.e., for each xo e [x n = 0} and T > 0, there is a constant L which is depending on 
Xq such that 

(5.16) \w(x, T) - Lx n \ < Co\x - x \ 1+y , in B^(x ) 
for some constant Co > 

Proof. For each T > 0, let r = r(lT) be given in Lemma 1531 for 1 = 1- (16)"~ 2 . We denote by 

Mo = max wq and mo = min wq 

K K 
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for the set K = [x € B+ : |x ; - P|,;| < f , Pj = (0, • • • , 0, \) e R" and i = 1, • • • , n). Then, by the Lemma |5jTJ 
there exists a solution / of 

1 / = on <9B+ 

with the conditions i£>(x, 0) > f(x, 0) in B+ and 

|V/(x, 01 > c > 0, W > ZT, x e <?Bi n {x H = 0}. 

16 

Hence, 

|Vw(x, 0| > |Vf(x, 01 > c > 0, Vf > ZT, x e dB\ n {x„ = 0) 

16 

and this implies that there exists a constant A > such that 

(5.17) w(x, > Ax„ in B+ x [ZT, T]. 

16 

In addition, there also exists a constant B > such that 

(5.18) w(x, < Bx„ in B^ X [ZT, T] 

16 

since Wq e C 0,1 (H) and the function Bx n is a solution of ( 15. It . Now we are going to show | |5.16l l with scaling 
properties. To get them, we define the new function by 

w'(x', f) = 16w(x, (x' = I6x, t' = (16) 2_It (t - IT)) . 
Then, the function w' is also a solution of j5.1|l with initial data w(; IT). By (|5.17t and l|5.18t , we also obtain 

(5.19) Ax' n < w'(x', t') < Bx' n , in B, + x [0, T]. 
Let's assume that Aq and Bo are the optimal constants for | |5.19| |. If 



x'eB,t :iv'(x',0)> 



Aq + Bq 



then there exist constants Si > and m.\ > such that 

(5.20) w'(x' r 0) - A x'„ > mi, in K. 

This is because w'(x', 0) — AqX„ is Lipschitz continuous with respect to space variable. On the other hand, if 



>- 2 W\, 



then there exist constants 6j > and m' Y > such that 



B x' n - w'{x', 0) > m[, in K. 

Note that the constants ni\ and are only obtained by \A - B \ and Lipschitz continuity We first assume 
that 



<A + B \ , 



>-\\n\ 



lx' eB r + :w'(x' r 0)>(- 
Then, by similar argument as above, there is a function f such that 

w'(x', t') - A x' n > f'(x', t') in B+ x (0, T) 

and 

\Vf'(x', t')\ > c > 0, W > IT, x' 6 <9B^ n \x' = 0}. 

16 

Note that m\ in l|5.20b is only depending on |B - A\ and Lipschitz continuity of if'. Hence, the constant Co is 
proportion to \Aq - Bq|. Thus, there exists a constant < k < || such that 



o/(x', f') - Ax,', > k(B - A )x,', in B+ x [ZT, T]. 



Similarly, if 
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then, we also have 

B x' n - w'{x', t') > k(B - A Q )x' n in B\ x [IT, T]. 

16 

Therefore, we get B 1 > A 1 > such that < A < A x < B 1 < B Q < co, B 1 - A 1 < (1 - k)(B - A ) and 
(5.21) Aix' n < iv'(x', t') < Bix' n , in B\_ x [IT, T]. 

16 

Hence, from 115.211 , 

A x x n < w{x, t) < B x x n , in B + r x [IT + (1 - 1)IT, T], 

By iteration arguments, we conclude that there exist series {A k } and {B/ c } such that < Aq < A\ < ■ ■ ■ < A k < 
A k+1 < ■ ■ ■< B k+1 <B k <---<B 1 <B <oo, B k+l - A k+l < (1 - x)(B k - A k ) and 

A k x n < w(x, t) < B k x n , in x [(1 - (1 - l) k+1 )T, T]. 

This will imply 

(522) t^l±L£5A< Bk -A k <(l-K) k (B Q -A ), inB^x[(l-(l-/)' c+1 )T,T] 

X n 16*+1 

for C = limjt-^ A k = lim/c-,,*, B k . Now we get J5.22II when x e B + , . Hence, by the the property of logarithm 
function, we get 

\w(x,T)-Cx„\ (16\x\\- l0S ^ (1 - K) 



inB+ 



Since < k < j|, we have < — log 16 (l — k) < 1. Hence 05.16b holds for constant y = — log(l - x) and the 
lemma follows. □ 

From the previous Lemma l5T2l we can cover the remaining part of missing step for the long time existence. 
Corollary 5.3. Under the hypotheses ofLemma \5.2\ we also have 

\w z (x) -L\< C\\x - x \ y , in B+(x ) 

for some constant Ci > 0. 

Proof. Let xq be a point on [x n = 0} and let X s = (x' r s) be a point in B+(xo)- Then, 

zv(x b ) — zv(x a ) 

(5.23) io z {x a ) = lim K — -. 

b^a b — a 

By the previous Lemma (Lemma l5.2l l, we can get 

(5.24) \w{x*) - La\ = C^x" - x | 1+) ' + h.o.t. of \x" - x \ 
and 

(5.25) \w(x b ) - Lb\ = d\x b - xo\ l+ y + h.o.t. of \x b - x \ 

for some constant C\ > if x" is sufficiently close to x b . Here, h.o.t. means 'higher order terms'. Then, by 
l|5T23l , H5.241 and l|525|, we obtain 



\zv z (x a ) -L\< 



lim W ^ ~ Lb ~ ~ ^ 



b^n b — a 

Since x" is chosen arbitrary in B^(xq), we can get a desired result. 



< d\x a - x \ y + h.o.t. of \x a - x \. 
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6. Proof of Main Theorem 



6.1. Local coordinate change. To motivate the proof of Main Theorem, we will first compute the transfor- 
mation of the equation 

ft = mf a Af 

when one exchanges dependent and independent variables near the boundary. This change of coordinates 
converts the boundary into a flat boundary. More precisely assume that the function f belong to the space 
C 2+a (U X [0, T]). Let X = (x°, f ) = (x'°, x° n , f ) at the boundary dU x (0, T]. We can assume (by rotating the 
coordinates) that f Xj (Xo) = ■ ■ ■ - f Xn i (Xo) = 0, and f Xrl (Xo) > 0. It follows from the Implicit Function Theorem 
that if the number r\ is sufficiently small, we can solve the equation z = f(x, t) with respect to x n around the 
point X . This yields to a map x n = h(x',z,t) defined in a small box B+(x'°, 0, to)- We wish to compute the 
evolution equation of h from the one of /. To compute the evolution of h, we use the identities 

fx„ = y z ' f x < = ~Th' ft ~ ~¥ z ' f XnX " = ~~k^ ZZ/ 
fxm = -jt z (j%h zz - 2^h XjZ + h x J (z = 1, • • • , n - 1). 

This follows immediately from above computations, the equation 

ft - mfAf = 

transforms into the equation 

1 + \V x >h\ 2S 



(6.1) 
Set 



M(h) = h f - z c 
_ 1 + |V T ,fr| 2 



a nk ' = 



2^ 

" h z 



and a kV = d^r 



for k', V = 1, ■ ■ ■ ,n — 1. Then, h is a solution of the equation 

L[w] = w t -z a (a kl w kl ) , (k, I = 1, ■ ■ ■ , n). 

The operator M defined above becomes degenerate when z = 0. We can also easily compute its first 
derivatives with respect to x„ (z = 1, ■ • ■ , n - 1): 



( 



(6.2) 
Set 



L[w : Xj] = w t - z a 



l + \V x ,h\ 2 \ v^ 2/l *; 



+ z a V x . -(V x .w) 



a 



1 + \V x ,h\ 2 



Jtfn 



h z 



a" k ' = and a k ' v = b m , 



for k', I' = 1, • ■ • , n - 1. Then, h Xj is a solution of the equation 

L[w : x { ] =w f - z a V, (a k %w) , (k, 1 = 1, • • ■ ,n). 
We finish with the proof of Theorem ll.il 

The Proo f o f Theorem \l.l\ By the Theorem 12. 11 there exists a solution / of £L3} which is C 2 '^ -smooth up to 
the boundary, when < t < % for some constant t > 0. Let (0, To) be the maximal interval of existence for 
Cg'^-solution / of ( |1.3t and suppose that To < °o. By scaling time around t = To, we can assume that 

T = l. 

We will show the solution / is in C 2 '~(D.) at t = 1. Then, by applying Theorem 12. 1 1 again, we can extend 
Cg'^-smoothness of / to f = 1 + Ti for some constant Tj > 0. By maximality of To = 1, contradiction arise. 
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Hence, To = °° and theorem follows. 

Since Q is compact domain, we can express Q as the finite union 



j \ 

(6.3) Q = Q U 



U>i 



of compact domains in such a way that 



dist(Q ,<?Q) > | > 



and for all I > 1 

Q; = B p (xi) n Q 

with B p (xi) denoting the ball centere at x\ e <9Q of radius p > 0. The number p > will be determined later. 
The equation 1 11.21 , when restricted on the interior domain Go, is nondegenerate. Therefore, the classical 
Schauder theory for linear parabolic equations implies that 

(6.4) feCf~(Q ) atf = l. 

Hence, we are going to concentrate our attention on the domains Q;, (I > 1), close to the boundary of Q. 
Under the local coordinate change above, we can change the function f defined on Q; into the function h 
defined on T)\ e H. In addition, 

h e C 2 /(Di) (0 < t < 1) . 

By the maximality of To = 1, h is no longer in Cg '' (£)/) at time t = 1. However the existence theory gave 
h e W 12 (D ; ) (See chapter 5 in for the details). By the Lemmas |3~T1 l3~2l and [331 the coefficients of L in gT) 
and J6.2I I satisfy the hypotheses in the corollary [53] and the Lemma [4.41 respectively. Then, the Lemma [4.41 
and the corollary 15. 31 tell us that 

Vh e C r s {Q+{x u 1)) where y - ^ and r < 1. 
Hence, the coefficients of the equation | |6.1| | satisfy the assumption of Lemma 3.2 in (6)- By Lemma 3.2 in J6), 

(6.5) /ieCf(Qtfe,l)). 

2 

Let's choose p > so small that £); e B\ (x{). Then 



heC 2 s ' r (Di) at r = l. 



This immediately implies 
(6.6) 

Therefore, by (16.311 , j6.4b and (|6.6b , we get 



(6.6) / e Cf r (£li) at t = 1. 



/ 6 C 2 s ' y (Q) at t = 1 



and theorem follows. 
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